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ABSTRACT 

Fluctuating velocity and temperature data collected from instru- 
ments mounted at three levels on the R/V ACAIJIA were analyzed. The 
structure of atmospheric turbulence was examined on the basis of 
spectral distributions of wind and temperature fluctuations over the 
ocean. The results of this study were utilized to determine the 
validity of present formulations and prediction techniques for the 
momentum flux and the temperature-structure parameter over water. 

For the stability conditions encountered, the momentum flux 
calculations support present expressions for the constant flux layer. 
Proper directional deviation from prediction of the variation of e 
with height was observed. When corrections for stability were applied, 
proper adjustment toward the predictions were noted. However, com- 
plete adjustment was not evident. This suggests other factors as 

2 

wind wave effects may be important, calculations were in agree- 

ment with present formulations and prediction techniques. 
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I. INTRODUCTION 



Interdisciplinary investigations are being conducted at the Naval 
Postgraduate School to correlate optical propagation and atmospheric 
turbulence data in the marine boundary layer. Such information is 
necessary for design of future naval systems. In these investigations, 
personnel from the Meteorology, Mechanical Engineering, Oceanography, 
and Physics departments have been obtaining shipboard observational 
data to evaluate optical propagation in the near-ocean environment. 

The purpose of this study is to examine the structure of atmos- 
pheric turbulence on the basis of spectral distributions of wind and 
temperature fluctuations. Turbulence similarity expressions have been 
tested primarily over land. Only a few over water results have been 
obtained. Observational results from this study will provide a basis 
on which to evaluate and improve existing techniques in predicting 
over water optical propagation. 

During January 1974, an observational experiment was made aboard 
the Naval Postgraduate School research vessel, R/V ACANIA shown in 
Figure 1. This experiment was conducted in Monterey Bay approximately 
two miles offshore. The experimental location provided the near-ocean 
environment for turbulence data collection as well as the platform for 
ship to shore optical propagation data collection. Meteorological 
measurements of wind and temperature were conducted at three levels. 

The scope of this study concerned the examination of the turbulent 
fluctuations of velocity and temperature in the marine environment. 
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Figure 1. The R/V ACANIA 
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The examination included an evaluation of the height variation in the 

constant flux layer of (1) e (the viscous dissipation of turbulent 

2 

kinetic energy) and (2) (the temperature structure parameter). 
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II. THEORETICAL CONSIDERATIONS 



A. THE TEMPERATURE- STRUCTURE PARAMETER 

An important parameter in describing optical propagation is the re- 

2 

fractive-index-structure parameter, C^ . According to Wyngaard, et al 
2 

(1971), C^ can be related to its counter part, the temperature-struc- 

2 

ture parameter, C^ through a semi-empirical theory. The temperature- 
structure parameter has the form 

= (79 X 10"^ P/T^)"^ . (1) 

2 

C^ can be estimated by several methods of measurement. Two of these 
methods will be discussed here. 

At a separation distance, r , of the order of inertial subrange 
scales, the temperature structure function in a locally isotropic field 
has the form 

[T(x) - T(x + r)]2 = ^2/3 ^2) 

2 

C^ can be found by making simultaneous measurements of temperature, 

at two positions separated by the distance r . A second method for 

2 

determination of C^ is by its relation to the one dimensional tem- 
perature spectrum which has the inertial subrange form 

4>^(K) = 8 (3) 

where 3 is a constant and K is the streamwise component of wave 
number. 
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Indirect means of expression in terms of other meteorological 

parameters is through turbulence similarity theory. Surface-layer- 
similarity theory predicts that when the statistics of the mean and 
turbulent flow fields are properly non-dimensionalized, they become 
universal functions of only a stability parameter. Wyngaard, et al 
(1971) chose to use the Richardson number, Ri, as the stability para- 
meter in doing calculations from overland data. The Richardson 
number is defined by 

Ri = 8/T_(d , e/dz) . 

(du/dz) 

where g = acceleration of gravity, T = mean temperature, d9/dz = 
potential temperature gradient, and du/dz = mean wind speed gradient. 
It was shown by Wyngaard that the temperature structure function could 
be expressed as a function of the Richardson number, potential tempera- 
ture gradient and height. 

(d0/dz)^ f (Ri) (5) 

The results of Wyngaard *s calculations are shown in Figure 2, verifying 

the similarity theory prediction for unobstructed overland data. 

2 

In order to fully understand the relation of to turbulence, 

one must consider e and , the viscous dissipation of turbulent 

kinetic energy and the viscous dissipation of temperature, through 
theoretical relations. Both the balance expression for turbulent 
kinetic energy and Kolmogorovas hypotheses must be examined. These 
theoretical considerations are reviewed in the following paragraphs. 
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Figure 2 



B. BALANCE EXPRESSION FOR TURBULENT KINETIC ENERGY 

The turbulent kinetic energy balance expression has the following 
general form 
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Neglecting horizontal eddy transports, energy flux in the vertical can 
be expressed as follows 
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The four terms on the right hand side of the equation represent pro- 
cesses which will cause a time variation of turbulent kinetic energy. 
These processes are 

i) Mechanical production of turbulent kinetic energy from the 
mean flow 

ii) Archimedes work term or the energy released due to the effects 
of buoyancy 

iii) Viscous dissipation of turbulent kinetic energy 

iv) Pressure work and flux divergence of turbulent kinetic energy 
due to transfer by turbulent motion 

Further simplification of the kinetic energy balance expression can be 

made by making several definitions and assumptions: 

i) Assuming that the turbulent kinetic energy is in a statisti- 
cally steady state, 

I? ^ 
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ii) Assuming horizontal homogeneity (i,e., the spatial variations 
of mean fluctuations are negligible), 

= 0 
3x 3y 



iii) Specifying the co-ordinate system such that 



“2 " “3 ' 0 



iv) Assuming the divergence terms are negligible 



r = n 

3x3 ^^3 ^ 2 ^ ® 



The turbulent kinetic energy balance expression can then be expressed 
in the following simplified form, for the case when the mechanical pro- 
duction of turbulent kinetic is equal to the viscous dissipation plus 
energy released or gained due to buoyance effects 

3u- 



—V V = C ^ - V 0 

13 3x^ ^ T 3 

3 o 



( 8 ) 



In cases where the atmosphere is in near-neutral stability conditions, 
the second term on the right hand side of the equation is negligible 
and Equation (8) can be written as 



2 ^ 
t 3x. 



= £ 



(9) 



where u^ = , The following expression for the vertical gradient 

of the mean wind also applies to the neutral condition 

3u^ 



3x. 



Kz 



( 10 ) 
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where K is a universal non-dimensional constant (i.e., the Von Kaman 
constant) and Z is height. Combining Equations (9) and (10) yields 
an expression to estimate u^ based on measurements of e at a known 
height Z 



C. VARIATION OF e WITH HEIGHT 

Equation (11) can be used to verify the constant flux assiamption. 
A straight line plot, with a -1 slope, of Ine against InZ implies 
u^ is a constant with height. For non-neutral stratification in the 
atmosphere, one must retain the buoyancy term in Equation (8) . Equa- 
tion (11) can be modified as follows to include the buoyancy effects. 



3 



u 



e 



KZ 



( 11 ) 



3 




( 12 ) 



or 



3 



^ KZ ^*^2 



(13) 



where 



L = -Tu^^/gK v^Q (Monin-Obuhkov Length) 

v^O = heat flux 

Ri = Richardson number 

= empirical function of Z/L 
= empirical function of Ri 



Equation (13) can be written in the following form 

3 



In e = -In Z + In (|) 2 (Ri) + In (~j^) 



(14) 
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From overland experiments, various expressions for (f)^ and (p^ have 
been derived. For cases of neutral stability, very little variation 
from a -1 slope is noted in a plot of In e versus In Z • In cases 
of instability, an effect on the plot of In e versus In Z is to 
increase the slope (in a negative sense), to greater than -1, Figure 3 
Illustrates these effects. 



D. VARIATION OF WITH HEIGHT 

2 

Equation (5) expressed as a function of the Richardson number, 

potential temperature gradient and height. The potential temperature 
gradient can be expressed in the following manner 



dZ ” Z '•L'' 



(15) 



where 



T* = v^e/Ku* 

L =* Monin-Obuhkov Length 

= empirical function of stability 



The Richardson number can also be expressed as a function of Z/L in 
Equation (5) 



Ri 



(16) 



Combining Equation (5) and Equation (16) results in the following ex- 

c n 2 

pression for 



_ 2 ^ 2 -2/3 Z 

Ct T* Z g^ (j^) 



(17) 
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Figure 3. An Illustration of the Effect of Stability on the 
Variation of e with Height 
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where is an empirical function of stability derived from overland 

data. For cases of neutral stability, = 1 , and 




(18) 



A plot of In versus In Z would yield a straight line with a -2/3 

slope • 

2 

Under stable conditions Wyngaard noted that decreased with 

-2/3 

height more slowly than Z . For unstable conditions, overland 

results indicate that g^ approximates the function 



g r A (l)-2/3 

83 3 



(19) 



Also, for unstable conditions. Equation (17) can be expressed as 




I T.2 z-4/3 



( 20 ) 



A plot of In versus In Z would yield a straight line plot with 

a slope of approximately -4/3, assuming overland results are valid 

2 

over water. Figure 4 shows the predicted variation of with 

height for various stability conditions. 



E. KOLMOGOROV’S SPECTRAL SIMILARITY EXPRESSIONS 

2 

The parameters e and can be estimated on the basis of 

Kolmogorov’s expressions for universal spectra (in the inertial sub- 
range) of both fluctuating velocity and temperature spectra. The 
inertial subrange in the velocity spectrum represents a range of eddy 
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Figure 4, 



An Illustration of the Effect of Stability on the Variation 
of with Height 
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sizes where energy is transferred from the production region to the 
dissipation region, (i.e., from lower to higher wave numbers). 

Figure 5 from Lumley and Panofsky (1964) shows the separate regions 
of a typical turbulence spectrum. 

According to Kolmogorov’s second hypothesis, the following func- 
tional relations hold 



♦ (K) - 


(21) 




(22) 



where ({)(K) and represent one dimensional velocity and tem- 
perature spectra and K represents the wave number. and are 

2 -1/3 

empirical constants. is defined to be e ^0 » therefore. 

Equation (22) is identical to Equation (3). Equations (21) and (22) 

2 

can be used to find estimates of e and by examining respective 

velocity and temperature spectra. 

Turbulent measurements are made at a fixed point in the flow and 
hence spectra are obtained for temporal frequency f . In order to 
use Equations (21) and (22), time and space scales must be related 
using Taylor’s (1938) ’’frozen turbulence” hypothesis 

K = 27Tf/u (23) 

where u is the mean wind speed. The hypothesis implies that turbu- 
lence remains unchanged during the time required for it to sweep 
past the probe. The following forms of Equations (21) and (22) feature 
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Figure 5. A Schematic Drawing Showing Spectral Transfer and 

Illustrating (A) Energy Containing Region (B) Inertial 
Subrange (C) Dissipation Range and (D) Universal 
Equilibrium Range 
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the transformation from frequency to wave number 



f (J)(f) = K (J)(K) = Z K ' 
f (J)^(f)= K (J)^(K)= C 2 



2/3 -2/3 



(24) 



(25) 



A second method for estimating e results from examining the 
upper limit of the inertial subrange of a velocity spectrum. In 
statistical equilibrium, a spectrum of turbulence is independent of 
the source of the turbulence and depends only on £ , Y ^ 

where Y is the kinematic viscosity* The following expression is 
applicable for the inner scale 




(26) 



where (lo) 



-1 



is proportional to the dissipation wave number. The 



dissipation wave number can be obtained from spectra which show scale 
separation between the inertial subrange and viscous dissipation 
regions* £ can then be determined from Equation (26). 
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III. OBSERVATIONS 



A. DATA COLLECTION 

Observational experiments were made from the R/V ACANIA anchored 
in Monterey Bay approximately two miles offshore. Observations were 
made at three levels on the forward mast and bow of the ACANIA. 

Figure 6 shows the mounting arrangements during the experiments. 

B. EQUIPMENT 

1. Fluctuating Velocity and Temperature Measurements 

Velocity sensors were (TSI) Model 1210 ”hot wire” probes with 
sliding support shields. Calibration of the sensors was made prior 
to the experiments. Band filtering was applied to the signal with cut- 
off frequencies of 2 hz and 2 khz respectively. The signals were then 
amplified and recorded on a 14 channel Sangamo Model 3562 FM recorder. 

Temperature sensors used were (TSI) Model 1210 platinum wire 
probes with sliding support shields. Electronic conditioning similar 
to that applied to the velocity signals was applied to the temperature 
signals. The temperature data was analog recorded simultaneously with 
the velocity data. 

2. Mean Velocity, Temperature and Humidity Measurements 

Mean velocity measurements were made at each level using C. W. 
Thornthwaite cup anemometers. The cup revolutions were counted by 
Hewlett Packard 5221A electronic counters. Direct readout of mean wind 
data were then made and recorded in a log. Mean temperature measure- 
ments at each level were made using Hewlett Packard quartz thermometers. 
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Figure 6. Mounting Arrangements 
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The signals were passed to a printer, allowing the temperature data to 
be printed out in a time series. Humidity data were collected in a 
fashion analogous to temperature data using a Hydrodynamics Digital II 
system and a printer. 

Figures 7 and 8 show the electronics involved in processing 
both mean and fluctuating data in the ACANIA^s laboratory. 

C. CALIBRATION 

1. Calibration of Hot Wire Sensors 

Velocity fluctuation measurements were made with hot wire sensors. 
These velocity sensors were calibrated before each experiment using a 
manometer and TSI 1125 Calibrator in accordance with procedures outlined 
by Johnston (1974). For any given mean wind speed, there is a corres- 
ponding voltage output. The relationship between wind speed and voltage 
is described by the following equation 

= a(u)^^^ + b (27) 

where v is voltage and u is the wind speed. Constants a and b 
are calibration curve slope and intercept respectively, and are determined 
experimentally . 

To analyze the recorded velocity data it was necessary to deter- 
mine a calibration factor such that 

c X v^ = u* (28) 
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Figure ?♦ Equipment Used Aboard R/V ACANIA to Obtain Fluctuating 
Velocity and Temperature Data. 
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Figure 8. 



Equipment Used Aboard R/V ACANIA to Obtain Mean Velocity 
and Temperature Data, 



31 





M 



where c = calibration factor, v* =» voltage fluctuation and u* = 
velocity fluctuation. Differentiation of Equation (27) yields 



r4vu 



1/2 



-] V* = u* 



( 29 ) 



where c , the calibration is then determined by the relation 

1/2 



[ 



4vu 



“] 



(30) 



For analysis, determination of the calibration factor consisted of 
(1) obtaining the slope of the calibration curve (2) selecting a mean 
wind speed value observed during the experiment and (3) entering the 
calibration curve with this mean wind value to determine a correspond- 
ing voltage. Figure 9 is an example of an actual hot wire calibration 
curve. 

2. Calibration of Platinum Wire Sensors 

Temperature fluctuation measurements were made using platinum 
wire sensors. A small temperature fluctuation causes a corresponding 
resistance change in the platinum wire. The resistance of a given wire 
can be related to temperature through the following equation 



R = Ro (1 + a (T-To)) 



(31) 



where R = resistance, T = temperature, To = 0®C, Ro = resistance 
at 0®C, and a = a constant for platinum. The resistance, in ohms, can 
be measured at any given temperature and thus Ro determined for each 
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Figure 9. Hot Wire Calibration Plot 
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MOT WIRE PROBE NUMBER 1 



wire. Differentiating Equation (31) results in 



dR = a Ro dT 



( 32 ) 



or 



R* = a Ro T* 



(33) 



which relates a small temperature change to a small change in the 
resistance of the wire. During analysis, resistance fluctuations are 
recorded as voltage fluctuations, through the electronics involved. 

A transfer function, which relates resistance to voltage, must be 
determined during calibration, i.e.. 



H , the transfer function can be found by substituting a variable 
resistor in place of the wire and recording voltage changes as a 
function of corresponding resistance changes. By plotting these points 
on a graph, the resultant slope will determine the transfer function 
H . Figure 10 is an example of an actual calibration curve used during 
one of the experiments. Substituting Equation (34) into Equation (33) 
results in an expression to relate recorded voltage fluctuations to 
temperature fluctuations 



V* = HR’ 



(34) 




(35) 



or 




(36) 



where C is the calibration factor. 
H 
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Figure 10, Platinum Wire Transfer Function Plot 
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IV. DATA REDUCTION AND ANALYSIS 



A. PRELIMINARY ANALYSIS 

Data period selection was based on evaluations of strip chart 
records of the signals. Criteria for selection of data periods con-- 
sisted of (1) signal strength (2) possible averaging time and (3) 
availability of mean profiles. Typical signals meeting the above 
criteria are shown in Figure 11. 

B. DIGITAL METHODS 

1. Analog to Digital Conversion 

Data from selected periods were prepared for further analysis 
by an analog to digital conversion of the magnetic tapes as described 
by Johnston (1974). The analog to digital process was performed on 
the hybrid computer system located on the fifth floor of Spanagel Hall 
at the Naval Postgraduate School. The hybrid system consists of an 
analog computer, COMCOR Ci 5000, interfaced electrically with a 
digital computer, XDS 9300. 

2. Seven Track to Nine Track Conversion 

The analog to digital process yields seven-track octal based 
records written on a digital tape. Analysis of the data was then per- 
formed on an IBM-360/67 computer. This computer requires nine-track 
hexadecimal based digital records. Conversions from seven-track to 
nine-track tapes were made with procedures described in Naval Post- 
graduate School Computer Center Technical Note #0211-08. Two programs 
used for this conversion were TAPEOUT and CONVERT. 
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Typical Velocity Signal 




Typical Temperature Signal 



Figure 11 
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3, Digital Analysis Programs 



Spectral analysis of the digital tapes was made using FTOR 
(Fast Fourier Transform) and SCOR (Spectral Analysis), programs which 
are available in the IBM 360/67 library. Descriptions of the FFT 
program package have been described by McKendrick (1972) and Johnston 
(1974). The output of th4 SCOR program consists of a tabular computer 
printout as shown in Figure 12 and a graphical plot of the spectrum 
versus frequency as shown in Figure 13, 

C . SCALING PROCEDURES 

During observational and analysis procedures, several gains and 
scaling factors were applied to the recorded signals. Prior to record- 
ing the data, an amplitude gain was applied to the signals. During 
the analog to digital process an additional gain was applied to the 
signal along with a scaling factor. The scaling factor resulted from 

the consideration that the maximum voltage input to the COMCOR Ci 5000 

23 

is 100 volts, and this corresponds to a digital maximum of 2 on the 
XDS 9300. Applying the respective gain and scaling factors resulted 
in a signal input to the IBM 360/67 as follows 

V (volts) X G1 X G2 X SC = v (volts IBM 360) (37) 

where v =* signal from velocity or temperature unit 

G1 = amplification prior to recording 
G2 = amplification during A-D process 
SC = analog to digital scaling factor 
The spectral analysis programs yielded plots of the energy density 
function <J)(f) , versus frequency f , The corresponding units 
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Figure 12. SCOR Spectral Statistics 
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Figure 13. SCOR Spectral Plot 



40 



2 

associated with c()(f) and f were VOLT SEC and 1/SEC respectively. 

Multiplying each point of the energy density function by its correspond- 

2 

ing frequency 4^(f)f results in a quantity with the units of VOLTS 
Taking the square root of this quantity, one can relate VOLTS to 
VOLTS’ (IBM 360) through Equation (37), Final scaling was accomplished 
by using Equations (29) and (35) which relates volts to either velocity 
or temperature fluctuations respectively. 
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V. INTERPRETATION OF SPECTRA 

A. SINGLE POINT METHOD 

The Single Point Method is perhaps the most straightforward method 

2 

used to obtain e and from velocity and temperature spectra. 

From the spectrum, a wave number, K , within the inertial subrange is 

2 

chosen. Using Equations (21) and (22), calculations of e and 
can then be made directly. 

B. INTERCEPT METHOD 

A second method described by Johnston (1974) can also be used to 

2 

evaluate e and from velocity and temperature spectra. This 

method has an advantage over the single point method in that many points 
are used to evaluate the spectral estimate. This is accomplished by 
taking the average slope through many points and determining an inter- 
cept value as discussed in the following paragraphs. 

1. Determination of e 

Equation (24) can be written in the following form 

In K(f)(K) = -2/3 In K + In (38) 

results in a curve with a slope 
At the point of intercept, 

(39) 



Plotting In K(j)(K) against -2/3 In K 

2/3 

of +1 and an intercept of In C^e 
In K equals zero and 



In = In C^e 



2/3 
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Letting A equal In K(f)(K) , the intercept value, e then becomes 



A 




3/2 



( 40 ) 



By determining the intercept value, e can be estimated* Figure 14 

shows a typical plot of this type using velocity data at three levels* 

2 

2, Determination of 

2 

A similar procedure can be used to determine . Equation 

(25) can be written as follows 

In K(J)^(K) = -2/3 In K + In C 2 0 ^ (41) 



Using the same procedure as described in the determination of e , 

2 

can be evaluated at the point of intercept in the following manner 



C 



T 



2 




(42) 



where B = In Kcf)^(K) at the ordinate intercept. Figure 15 is an 
example of an intercept plot using temperature data at two levels* 



C* INNER SCALE METHOD 

A third method of estimating e was to consider the inner scale. 
Equation (26) describes an expression relating e to the separation 
between the inertial subrange and dissipation region of the spectra. 

By choosing the wave number at which the dissipation region occurs in 
the velocity spectra, e can be evaluated* This method was quite sub- 
jective with regard to determining the exact point of scale separation* 
However, the results were qualitatively good when compared to the 
intercept method. 
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Figure 14. Intercept Plot for Velocity Data at Three Levels. 
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Figure 15. Intercept Plot for Temperature Data at Two Levels. 
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VI. RESULTS 



A. SUMMARY OF PERIODS EXAMINED 

Fluctuating velocity and temperature data were collected from three 
levels aboard the R/V ACANIA on 18 January 1974, Selected periods for 
analysis were chosen for three time intervals between 1834 and 1932, 

Table I lists periods examined, heights at which data were collected, 
parameters measured and stability conditions for each period. 

B. VELOCITY RESULTS 

Velocity spectra were obtained for simultaneous measurements at three 
levels during each period. From these spectra e and u^ were computed 
by methods described in Sections V-B-1 and V-C. Table II lists e and 
u^ calculations for each of these periods. 

The original experimental plans were to measure mean wind data at 
each of the levels. Due to problems encountered with the read-out counter, 
the wind data appeared to be in error. The problem was resolved by use of 
mean wind measurements at one level, the ACANIA^ s aerovane, and then esti- 
mating a vertical wind field based on the Monin-Obuhkov prediction. A 
Fortran computer program, PROPAR, described by Cardone (1969) was used 
to estimate the vertical wind field. This program uses inputs of (1) 
mean wind speed (2) anemometer height (3) air-sea temperature difference 
and (4) thermometer height to compute the output parameters. These com- 
puted parameters included (1) u^ = friction velocity (2) Zo = roughness 
parameter and (3) L = Monin-Obuhkov Length. Mean wind was calculated 
at each level using the following equation 
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TABLE I 



Data Periods Considered in This Study 
18 January 1974 





Time 


Height 


Parameter 


Stability 


(local) 


(cm) 




(Ri) 


1834-1854 


394.0 
(level 1) 


U,T,Q,U',T* 


-0.136 




658.0 
(level 2) 


U,T,Q,U',AT' 






822.0 
(level 3) 


U,T,Q,U',T* 




1855-1913 


394.0 
(level 1) 


U,T,Q,U',T’ 


-0.293 




658.0 
(level 2) 


U,T,Q,U',AT' 






822.0 
(level 3) 


U,T,Q,U',T* 




1914-1932 


394.0 
(level 1) 


U,T,Q,U',T’ • 


-0.923 




658.0 
(level 2) 


U,T,Q,U',AT* 






822.0 
(level' 3) 


U,T,Q,U',T' 
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TABLE II 



e and Momentum Flux Calculations for 18 January 1974 using 
the Zero Intercept Method 



Height 



t 2 -3, 

(cm sec ) 



/ “In 
(cm sec ) 



Time 1834-1854 



394.0 


83.43 


23.60 


658.0 


54.74 


24.33 


822.0 


22.83 


19.57 




Time 1855-1913 




394.0 


342.96 


37.80 


658.0 


79.21 


27.51 


822.0 


36.72 


22.93 




Time 1914-1932 




394.0 


112.28 


26.05 


658.0 


13.67 


15.32 


822.0 


22.22 


19.40 



48 



u 



(43) 



(i- <fe> - *1 <!>) 

where is an empirical function of Z/L described by Paulson (1970) • 

The purpose in defining the wind field distribution was to obtain advec- 
tion speeds for relating temporal and spatial frequencies. The need for 
wave number (spatial frequency) description was identified in Section 
II-E. 

Figures 16, 17 and 18 show typical velocity spectra at each level for 
one of the periods. In each case a -5/3 slope is observed for the 
inertial subrange region. At higher frequencies, the spectra exhibit 
greater slopes. The regions of different slopes represent different 
scales in the turbulent regime. The scale represented by greater slope 
at higher frequencies, is the dissipation region. Closer investigation 
reveals a decrease in the frequency where the slope changes as height 
increases. However, as mentioned earlier and as can be seen here, the 
exact point of scale separation between the inertial subrange and dissi- 
pation region is difficult to determine. 

Results on e versus height for the velocity data are shown in 
Figures 19, 20 and 21 where In e is plotted against In Z for each 
data period. Negative slopes, greater than -1, are noted in each case. 

As stated in Section II-C, this agrees with present boundary layer theory 
for unstable atmospheric conditions. 

The deviation of the slope from -1 was examined to see if it could 
be accounted for by the influence of the unstable stratification. The 

plots were corrected for stability using Equation (13) where (f >2 ~ 

-1/4 

(1 - 12Ri) , A slight shift toward a -1 slope was observed; however, 
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Figure 16. Velocity Spectrum Level 1 for 1855-1913, 18 Jan 74. 
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Figure 17. Velocity Spectrum Level 2 for 1855-1913, 18 Jan 74. 
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Figure 18. Velocity Spectrum Level 3 for 1855-1913, 18 Jan 74. 
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Figure 19. Measured Values of e vs Height for 1834-1854 
18 Jan 74. * 
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Figure 20. Measured Values of e vs. Height for 1855-1913, 18 Jan 74. 
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Figure 21. Measured Values of e vs. Height for 1914-1932, 18 Jan 74. 



complete adjustment to a -1 slope was not evident. Similar results 
were found by Johnston (1974), He states **The fact that the stability 
corrections did not completely restore the -1 slope, even in the 
stable case, leads to the examination of the possible effects of wind- 
wave coupling,” 

The above results imply that overland prediction for the influence 
of stability are not completely valid over water and that wind-wave 
coupling or wave related velocity fluctuations might be the reason. 
There is still not enough experimental and theoretical knowledge to 
describe the waves* influence upon the data. 

Several consistent features can be observed in the results. During 
the first two periods there appears to be a bending or ”kink” in the 
plots such that they are concave downward. In the third data period, 
the ”kink” is concave upward. However, strip chart records show a 
deterioration of the velocity signal at level 2 toward the end of the 
period. This apparently resulted in low estimates of e for that 
level. It is not felt that this ”kink” is due to the ship’s influence, 
Welsh (1974) found that ship related motion had little effect on turbu- 
lence measurements at these frequencies, Johnston (1974) did not find 
this characteristic bending in similar plots of In e versus In Z 
for unstable conditions. In these results, both the intercept method 
and inner scale method exhibit the ”kink”. Based on only two periods 
of analysis, conclusions as to the cause of this feature cannot be 
asserted, 

C, TEMPERATURE RESULTS 

Fluctuating temperature data were also obtained at three levels for 

2 

each period. At levels 1 and 3, was determined from the 
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one-dimensional temperature spectrum as described by Equations (3) and 

(22), At the mid level, two probes were separated by r = 10 cm and 

2 2 

determined by using Equation (2), Table III lists the calcu- 

lations for each of the periods. Temperature spectra at three levels 
are shown for one of the periods in Figures 22, 23 and 24, These are 

typical of the spectra for the other periods examined, 

2 

Results of versus height are shown in Figures 25, 26 and 27, 

2 

where In is plotted against In Z , The general slopes of these 

plots agree reasonably well with the -4/3 prediction for unstable con- 
ditions discussed in Section II-D, Periods 1 and 2 both show close 
agreement with the -4/3 prediction. However, period 3 shows a much 
' greater deviation from -4/3, Comparing Richardson numbers from Table I, 
it is noted that period 3 was much more unstable than the first two 
periods. One might expect slightly greater negative slopes for greater 

instability; however, -4/3 should be a limiting slope if overland results 

2 

are valid over water. Since is a function of both e and , 

the possible wind wave effects on e discussed in the previous section, 

2 

would also influence , This is a possible explanation for any 

2 

deviation in from a -4/3 slope, 

2 

Closer examination of In C^ versus In Z reveals "kinking** or 

bending similar to that discussed in the velocity results section. In 

both periods 1 and 2, this bending phenomenon is concave downward as it 

was for e calculations. During period 3 the general appearance of the 
2 

kink in the C^ plot approximates the appearance of the kink in the e 
plot, disregarding any influence of a deterioration in the velocity 
signal during that period. In this case however, the bending is in 
the opposite direction. 
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TABLE III 



C 

T 



2 



Calculations for 18 Jan 74 using the Zero 
Intercept Method 



Height 

(cm) 






Time 

394.0 

658.0 

822.0 



1834-1854 

.00312 

.00239 

.00168 



Time 

394.0 

658.0 

822.0 



1855-1913 

.00600 

.00558 

.00428 



Time 

394.0 

658.0 

822.0 



1914-1932 

.00934 

.00416 

.00393 
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Figure 22. Temperature Spectrum Level 1 for 1855-1913, 18 Jan 74. 
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Figure 23. Differential Temperature Spectrum for Fluctuating 
Temperature Data Separated by r = 10 cm, Level 2 
for 1855-1913, 18 Jan 74. 
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Figure 24. Temperature Spectrum Level 3 for 1855-1913, 18 Jan 74. 
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Figure 25. 



2 

Measured Values of C„ vs. Height for 1834-185A 
18 Jan 74. ^ 
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2 

Figure 26. Measured Values of C„ vs. Height for 1855-1913, 
18 Jan 74. 
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Figure 27. 



2 

Measured Values of C vs. Height for 191A-1932 
18 Jan 74. 
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The above results give further support to the consideration that 
the bending is actual and not a result of random errors within the 
system, since these spectra were obtained from two independent sources, 
velocity and temperature measurements. However, no further explana- 
tion as to the cause of the bending is offered due to the limited 
amount of data considered. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



Observed height variations of e and C^ tend to support present 

formulations and prediction techniques. However, several exceptions 

were noted and will be summarized in the following paragraphs. 

For unstable atmospheric conditions slopes of In e versus In Z 

were greater than -1. When these slopes were corrected for stability, 

there was a tendency for the slopes to shift towards a -1 slope. The 

observation was made that complete adjustment to a -1 slope was not 

evident, suggesting other factors such as wind wave effects may be 

important in near-ocean boundary theory. The observation was also made 

that In e versus In Z was not a linear plot. There appeared to be 

bending or ’’kinking’* which is not predicted in present formulations. 

No definite reason was given to account for the apparent kinking due 

to the limited amount of data examined. 

2 

A plot of In C^ versus In Z resulted in -4/3 slopes for unstable 

conditions. Some slight deviations were noted and attributed to other 

factors mentioned in the previous paragraph. The characteristic ’’kinking” 

2 

or bending in the In C^ versus In Z plots was also evident. 

Spectral analysis appears to be a satisfactory method for determining 

2 

e , momentum flux and C^ estimates over the ocean. However, digital 
methods used are very time consuming due to the amount of tape handling 
and large computer programs used in spectral analysis. A real-time 
spectrum analyzer would increase the amount of data that one could 
examine due to the time saved in the analysis. It is suggested that 
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subsequent fluctuating velocity and temperature experiments aboard the 
ACANIA be examined using a spectrum analyzer in the analysis process. 
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